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We theoretically study two atoms with p-wave interaction in a one-dimensional waveguide, and 
investigate how the transverse anisotropy of the confinement affects the two-body state, especially, 
the properties of the resonance. For bound-state solution, we find there are totally three two-body 
bound states due to the richness of the orbital magnetic quantum number of p -wave interaction, while 
only one bound state is supported by s-wave interaction. Two of them become nondegenerate due 
to the breaking of the rotation symmetry under transversely anisotropic confinement. For scattering 
solution, the effective one-dimensional scattering amplitude and scattering length are derived. We 
find the position of the p -wave confinement-induced resonance shifts apparently as the transverse 
anisotropy increases. In addition, a two-channel mechanism for confinement-induced resonance in 
a one-dimensional waveguide is generalized to p-wave interaction, which was proposed only for s- 
wave interaction before. All our calculations are based on the parameterization of the 40 K atom 
experiments, and can be confirmed in future experiments. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.50.-s, 67.85.-d 


I. INTRODUCTION 

Nowadays, the confinement-induced resonance (CIR) 
as one of the most intriguing phenomena in low¬ 
dimensional systems has attracted a great deal of inter¬ 
est. It was first predicted by Olshanii when considering 
the two-body s-wave scattering problem in quasi-one- 
dimensional (quasi-ID) waveguides [1], Subsequently, 
this study was extended to quasi-two-dimensional (quasi- 
2D) systems [2, 3]. In the past decade, an impressive 
amount of experimental and theoretical efforts have been 
devoted to confirm the existence of CIRs and explore the 
important consequences [4-18]. To date, CIR has already 
become a fundamental technique in studying strongly in¬ 
teracting low-dimensional quantum gases. 

The p -wave interaction is of particular interest, since 
it is the simplest high-partial-wave interaction with non¬ 
zero orbital angular momentum 1 = 1. This leads to a 
spatial anisotropic scattering in few-body physics, and 
open a new method to manipulate resonant cold atomic 
interactions by using the magnetic field vector, as indi¬ 
cated in our previous work [18]. For many-body physics, 
the high-partial-wave interaction may also result in phys¬ 
ically abundant quantum phase transitions, which are 
absent in s-wave interaction [19-22], In addition, the 
p -wave scattering dominates the low-energy interaction 
between spin-polarized Fermi atoms (all in one hyperfine 
state) due to the Fermi-Dirac statistics, which provide an 
ideal candidate to study the p -wave scattering properties 
in cold atoms [23, 24]. 
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For low-dimensional quantum systems, the external 
trap potential is one another freedom to manipulate two- 
body resonant interactions, and it is interesting to iden¬ 
tify how CIR is affected when the external confinement 
changes. For example, in ID Bose 133 Cs atoms with s- 
wave interaction, it is found that the transverse harmonic 
anisotropy shifts the position of CIR [12, 13], besides, 
CIR splits due to the anharmonicity of the transverse 
trap [10, 14, 15]. Naturally, it is also important to in¬ 
vestigate how the external confinement affects the low¬ 
dimensional resonant p -wave interaction between spin- 
polarized Fermi atoms, such as 40 K atoms. 

In this paper, we study the influence of the transverse 
anisotropy on two-body state with p -wave interaction, es¬ 
pecially, the scattering resonance in quasi-ID waveguides. 
We solve the two-body problem in a quasi-ID waveguide 
under transversely anisotropic confinement. The p-wave 
interaction is modeled by the contact pseudopotential in¬ 
cluding the effective range [25]. This model neglects the 
spatial anisotropy of p -wave interaction , which is not 
the focus of this work. However, it is already sufficient 
enough to study how the external trap affects the p -wave 
CIR. Due to the non-zero orbital angular momentum, 
we find there are totally three bound states for p -wave 
interaction in the waveguide, while the s-wave pseudopo¬ 
tential supports only one two-body bound state. Two 
of them become nondegenerate as the transverse con¬ 
finement changes from isotropic to anisotropic, which 
breaks the rotation symmetry around the waveguide. For 
scattering solution, we calculate the effective ID scat¬ 
tering amplitude as well as the ID scattering length, 
and predict the scattering resonance for any transverse 
anisotropy, whose position shows an apparent shift with 
increasing anisotropy. In addition, we present an effec¬ 
tive two-channel mechanism for p -wave CIR, which was 
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first proposed only for s-wave interaction [4] . In this two- 
channel picture, the transverse ground state and the re¬ 
maining transverse excited modes play roles of open and 
closed channels, respectively. When a bound state in the 
closed channel exists, and becomes energetically degen¬ 
erate with the scattering threshold of the open channel, 
a scattering resonance is expected. All our calculations 
are based on the parameterization of the 40 K atom ex¬ 
periments [7, 24], and can be confirmed in future exper¬ 
iments. 

In the following, we first study the properties of two- 
body bound states with p -wave interaction in a trans¬ 
versely anisotropic waveguide (Sec. II), and then present 
the scattering solution in Sec. Ill as well. We calculate 
the effective ID scattering amplitude and ID scattering 
length, and discuss how the transverse anisotropy of the 
confinement affects the p -wave CIR in details. In Sec. IV, 
a two-channel mechanism is presented for p -wave CIR. 
The main results are concluded in Sec. V. 


II. BOUND-STATE SOLUTIONS 


In order to simplify the problem, let us consider two 
atoms confined in a quasi-ID waveguide with a tight 
transverse harmonic potential (in x—y plane). The trans¬ 
verse trap frequencies are 0J x , y , and the atoms can freely 
move along z axis. The transverse anisotropy is char¬ 
acterized by p = LOx/ojy . Unlike the s-wave interaction, 
the p -wave interaction is strongly energy-dependence due 
to the narrow-width property [24]. Thus, one additional 
parameter named effective range should be included [25— 
28]. In harmonic traps, the center-of-mass (c.m.) motion 
is decoupled from the relative part, and we only need to 
solve the relative motion, while the c.m. motion is a sim¬ 
ple harmonic oscillator. By dropping the c.m. motion 
off, the relative motion of two atoms is described by the 
following Hamiltonian, 

fV d 2 

* = + ^ + (1) 

where 


Hi 


H 2 / d 2 d 2 \ 

2/i \<9x 2 dy 2 ) 


( rfx 2 + y 2 ) , (2) 


/i is the reduced mass, and w is the trap frequency in y 
aixs (we omit the subscript y without ambiguity). The 
interatomic interaction Vi (r) is modeled by the p -wave 
pseudopotential [25], 

Vl (r) = T (si “ U* 2 ) ^ s(r) ^ r ^’ (3) 


atoms as E = h 2 k 2 / (2p) . The symbol V (\^) denotes 
the gradient operator that acts to the left (right) of the 
pseudopotential. 

For the bound-state problem, the wavefunction can for¬ 
mally be written as 

ip (r) = - J d 3 r'G E (r, r') Vi (r') ip (r'), (4) 

where Ge (r, r') is the single-particle Green’s function 
with energy E , and satisfies 

+'H±-E S jG E {r,r') = 6(r-r'). (5) 

The Green’s function can be expanded in series of the 
eigen-states of non-interacting Hamiltonian as, 


OO OO n oo 

Ge (r, r ') = X X / dkz 

711=0 n 2 =0‘'~°° 

C nin2 {k z ) <p ni ) <Pn 2 (|) -j=, (6) 

where cp v (•) is the eigen-state of ID harmonic oscillator, 
d = sJhJyJjj is the harmonic length in y aixs, and k z is 
the wavenumber along the waveguide (z axis). Insert¬ 
ing Eq. (6) into Eq. (5) and using the completeness of the 
eigenstates of the non-interacting Hamiltonian, we may 
easily obtain the coefficients C nin2 ( k z ) , and then after 
some straightforward algebra (the derivation is similar to 
Eq.(24) in [12]), the Green’s function takes the following 
integral representation, 
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which is valid ior e = E/huj — (77 -1-1) /2 < 0 . Combining 
Eqs.(3) and (4), the bound-state wavefunction takes the 
form 


7 xh 2 f 1 1,2 

ip(r) = -- -rik 

/i \ui 2 


-1 






( 8 ) 


J r'=0 


where v\ and r\ are the three-dimensional (3D) p-wave 
scattering volume and effective range (with a dimension 
of inverse length), respectively, which can be tuned by 
using p -wave Feshbach resonances [23, 24], k is the rel¬ 
ative wavenumber, related to the relative energy of two 


where the summation is over n = x,y, z , and r x = 
x, r y = y, r z = z . Here, we have defined the coefhcent 






r—0 


( 9 ) 
























3 


Acting ■§pjf 3 -^ r - on both sides of Eq.(8) and setting r = 0 For any nonzero vector r , it is easy to show that 
, we obtain the following secular equation, 


E 


where 


7T H 2 / 1 1 2 \ 1 A r 

-- ~ r lK Ain ~ Oln 

n n = o, 
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,(13) 

(14) 

(15) 


Substituting Eq.(12) into Eq.(ll), it directly yields 

A/ti = ^ 3 / 2 ^ 5 y v j Fn ( e > °) Sln ’ ( 16 ) 

where Si n is Kronecker delta function. Combining 
Eqs.(10) and (16), the vanish of the determinant of the 
coefficient matrix in Eq.(10) determines the binding en¬ 
ergy, which yields 


quantum number as superposition of ip x (r) , ifj y (r) , and 
z (r) , i.e., 


V>m=±l (r) 
V>m=o (r) 



[ip x (r) ± ii/j y (r)] , 



(19) 

( 20 ) 


— = — j=Fi (e,0) +dn • E {I = x,y,z). (17) 

V\ \J 7T TUjJ 

Obviously, there are totally three bound states for p-wave 
interaction, while the s-wave pseudopotential supports 
only one two-body bound state. This is resulted from the 
richness of the orbital magnetic quantum number of the 
p-wave interaction. The corresponding (un-normalized) 
bound-state wavefunction is 


ipi (r) = n 


1 

^3 


E 1 
Huj rd 2 



■?i (e,r) 


(l = x,y,z). 


(18) 

Obviously, these bound states can be classified as two 
kinds by the transverse parity (in x or y axis), i.e., 
ij} x ,y (r) with odd transverse parity and ip z (r) with even 
transverse parity. For a transversely isotropic confine¬ 
ment (p = 1), there is a rotation symmetry around the 
z axis, and ip x (r) and ip y (r) are degenerate. Thus, at 
small separation, i.e., r ps 0 , we can construct another 
set of eigen-wavefunctions with specific orbital magnetic 


or explicitly, 
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( 22 ) 

with specific orbital magnetic quantum number m = 
0, ±1 , and 


Ti ± i (f) 


Yio (r) 


3 x ±iy 



(23) 

(24) 


are spherical harmonic functions. As the transverse con¬ 
finement becomes anisotropic, the rotation symmetry 
around the z axis is broken. This results in the non¬ 
degeneracy of ip x (r) and ip y (r) , and they are no longer 






















































4 



B{ G) 


Figure 1: (Color online) The bound-state energy of two 40 K 
atoms in the |.F = |,rriF = — ■§) hyperline state confined in 
a transversely isotropic waveguide near the p- wave Feshbach 
resonance centered at Bo = 198.8G in three dimension [7, 
24]. The dotted horizontal lines are the thresholds that the 
bound states merge into the continuum as the magnetic field 
increases. The magnetic field B dependence of the scattering 
volume ui as well as the effective range ri for 40 K is from 
Ref. [24], The insert shows the details of the curves in a smaller 
magnetic field range for eye clarity. 


the superposition of those with specific orbital magnetic 
quantum number at small r . 

We predict the binding energy of two 40 K atoms in 
the |F = | ,rriF = — |) hyperfine state as functions of the 
magnetic field near the p -wave Feshbach resonance cen¬ 
tered at B 0 = 198.8G in three dimension (3D) [7, 24], con¬ 
fined in a transversely isotropic waveguide (Fig.l) as well 
as in a transversely anisotropic waveguide (Fig. 2). For 
the isotropic confinement (p = 1), ip x (r) and ip y (r) are 
degenerate as we anticipate (see Fig.l), and they become 
distinguishable in energy due to the rotation symmetry 
breaking around the z axis when the transverse trap be¬ 
comes anisotropic (see Fig. 2). In addition, as the mag¬ 
netic field increases, these bound states gradually merge 
into the continuum at different energy thresholds, i.e., 
(377 + 1) /2 , (77 + 3) /2 , (r} + 1) /2 , respectively. 


III. SCATTERING SOLUTIONS AND 
CONFINEMENT-INDUCED RESOANNCES 
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Figure 2: (Color online) The same as Fig.l but for a trans¬ 
versely anisotropic confinement, i.e., p = 2 . 


ip(r) = ip 0 (r) - J d 3 r'G E (r, r') Vi (r') ip (r') , (25) 

where ipo ( r ) is the incident wavefunction. Since the en¬ 
ergy of the relative motion is just above the transverse 
zero-point energy, the atoms should enter from the trans¬ 
verse ground state, and the incident wavefunction takes 
the form of 

V’o (r) = ^0 (|) *sin(fcz 2 ). (26) 

Here, we have considered the exchanging anti-symmetry 
of two identical fermions. Remind that </>o (•) is the ID 
harmonic ground-state wavefunction and k z = \/2e/d . 
The two-body Green’s function G E (r, r') is the simple 
analytical continuation of Eq.(7) from e < 0 to e ~ 0 + . 
Substituting the p -wave pseudopotential (3) into Eq.(25), 
and at large separation, i.e., r —00 , we find the scat¬ 
tering wavefunction (25) behaves as 

■ ■ n \ ( \ ( 1 1 , 2 \ 1 vW /4 ^ ik PI 

1 sin (k z z) — sgn ( 2 )--qr —- 7 -e * l ~ l , 

\r;i 2 / d 

(27) 


In this section, let us consider the low-energy p-wave 
scattering problem in the ID waveguide with energy just 
above the transverse zero-point energy, i.e., e —> 0 + . 
According to the Lippman-Schwinger equation, the scat¬ 
tering solution of the Hamiltonian (1) can formally be 
written as 


where 1Z Z is defined as in Eq.(9). For this ID scattering 
problem in the waveguide, we anticipate the scattering 
wavefunction at large distance takes the form 


ip (r) - cpo 



0 ) i sin (k z z) + sgn (z) f p e lkz ^ z ^ 


( 28 ) 
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where f p is the effective p -wave ID scattering amplitude. 
Comparing Eqs.(27) and (28), we easily obtain the effec¬ 
tive ID scattering amplitude 


fp ~ 



1 , 1 vW /4 ^ 

2 r ' k ) - d - 


(29) 


The coefficient 7 Z~ is determined by the asymptotic be¬ 
havior of the scattering wavefunction (25) at small dis¬ 
tance, i.e., r ss 0 , and we find 
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2yf2 r 5 / 2 
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For the low-energy ID scattering problem, i.e., k z d <C 
1 , it is convenient to use the effective ID scattering 
lenght lip to characterize scattering resonances, which 
can be defined from the scattering amplitude at zero en¬ 
ergy (e « 0 +) as [ 8 ] 


fp — 


ik z 


ik z + l 


-l 1 
ID 


and 


h d _ R /- 
— = 6 ^ 


£-*,.2+1 +_Lr,( 0,0) 

Wl 2 V7T 


(32) 


• (33) 


The divergence of the ID scattering length character¬ 
izes the scattering resonance in the waveguide. By tuning 
the 3D scattering volume v± and the effective range ri 
according to the p -wave Feshbach resonance, the p -wave 
CIR can be reached in the ID waveguide when 



V! 


r]+l 

2 


- jr/( 0 , 0 ). 

V k 


(34) 


For a transversely isotropic confinement, i.e., r) = 1 , the 
ID resonance condition (34) returns to the well-known 
result [5, 8 ] 



Vi 



(35) 


where £ (•) is the Riemann Zeta function, and note that 
we include the 3D p -wave effective range in our expres¬ 
sion. 



Figure 3: (Color online) The p-wave ID scattering length 1\d 
as functions of the magnetic field B in the ID waveguide with 
different transverse anisotropy rj . Here, the p-wave interac¬ 
tion of two 40 K atoms in the F = | ,m,F = — |) hyperfine 
state is tuned by using the p-wave Feshbach resonance cen¬ 
tered at Bq = 198.8G in three dimension [7, 24], 


In Fig. 3, we present the ID scattering length lm 
as function of the magnetic field strength B in the 
waveguide with three typical transverse anisotropies rj = 
1, 5, and 10. Here, we still consider two 40 K in the 
| i 7 '=|,mp = —-() hyperfine state near the p- wave Fes¬ 
hbach resonance centered at Bq = 198.8G in 3D [7, 24] . 
We find that, as the transverse anisotropy 77 increases, the 
resonance position of p -wave CIR shifts to a higher mag¬ 
netic field strength. More apparently, we show how the 
resonance position denoted by the magnetic field strength 
B shifts with the transverse anisotropy 77 in Fig. 4. 


IV. THE EFFECTIVE TWO-CHANNEL 
MECHANISM 

For s-wave CIR in a quais-lD waveguide, there is a 
simple two-channel picutre first introduced by Bergeman 
el al. [4], and then extended to the case of the trans¬ 
versely anisotropic confinement [12]. We may understand 
this two-channel mechanism as follows: due to the low 
temperature and tight transverse confinement, the two 
atoms can only enter from the transverse ground state, 
while the transverse excited modes are all closed. There¬ 
fore, the transverse ground state and the manifold of the 
remaining transverse excited states may be regarded as 
the open and closed channels, respectively. If a molecu¬ 
lar state exists in the closed channel, a zero-energy scat¬ 
tering resonance occurs when this molecule energetically 
coincides with the continuum threshold of the open chan¬ 
nel. In this section, we aim to study whether this simple 
two-channel picture is still valid for p -wave interaction. 

Owing to the separability of c.m. motion and rela- 
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Figure 4: (Color online) The position of the p-wave 
confinement-induced resonance in the ID waveguide denoted 
by the magnetic field strength B < ' R ' 1 as a function of the trans¬ 
verse anisotropy g . In this plot, we tune the transversely 
anisotropy g by increasing uj x with fixed uj y . The correspond¬ 
ing parameters are still from the 40 K experiments [7, 24]. 


five motion in harmonic traps, we still only focus on the 
relative-motion Hamiltonian (1). According to the two- 


channel mechanism proposed for s-wave interaction [4], 
the total relative-motion Hamiltonian T~L may formally 
be splitted into three terms, i.e., "H op , T~L c i and W corre¬ 
sponding to the open channel, closed channel, and cou¬ 
pling part, respectively, 

n = Uop + Hd+W 

= P g 'HP g + P e UP e + (P g 'HP e + H.c), (36) 

where P g = |g) (g\ , P e = l a ) ( a l are the corre¬ 

sponding projection operators, and | g) , |a) are the trans¬ 
verse ground and excited states, respectively. In the fol¬ 
lows, we are going to show that the crossing of the molec¬ 
ular state in the closed channel T~i c i and the energy con¬ 
tinuum threshold of the open channel, i.e., Huj {g + 1) /2 
, denotes the position where the p- wave CIR occurs. 

In order to obtain the molecular state of the closed 
channel, we need to solve the Schordinger equation 
Pciipci ( r ) = -E'V’cZ (r) . The procedure is similar to that 
of obtaining the bound states of the full Hamiltonian (1). 
However, we need to project out the transverse ground 
state, and this means the transverse ground state, i.e., 
the n\ = TI 2 — 0 term, should be excluded in the expan¬ 
sion of the Green’s function, i. e., Eq.(6) . Similarly, we 
obtain the Green’s function for the closed channel Hamil¬ 
tonian, 


G * )(r - r ' ) = ?s 7 %k;“ p 
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-1 


(37) 


Then the bound-state wavefunction ip c i (r) takes the form 

-l 
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(38) 


r'—O 


where 


After straightforward and similar algebra as that in 
Sec. II, we find there are also three bound states in the 
closed channel, and two of them are with odd transverse 
parity, while that of the third is even, which takes the 
form, 


ipd,z (r) = 2 
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K (e>r) = 
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■ ( 41 ) 


The corresponding binding energy of the bound state 
with even transverse parity satisfies, 

d 3 6 . . £ 

— = — f K e,0) + dn-—. (42) 

ui V 7r nw 

Obviously, T' z (e, 0) takes the form of Eq.(31). If two 
atoms enter from the transverse ground state (open chan¬ 
nel) with energy below the transverse excited states 
(closed channel), they may only be coupled to the molec¬ 
ular state ip c i,z (r) with even transverse parity, due to the 
parity conservation. When this molecular state energeti¬ 
cally crosses the continuum threshold of the open channel 
(e ~ 0), a zero-energy scattering resonance occurs. From 
Eq.(40), the resonance condition is given by 



vi 


r/ + 1 
2 




J-'(0,0), 


(43) 


which returns to Eq.(34) again. The binding energy of 
the 40 K- 40 K molecular state 4> c i,z (r) in the closed chan¬ 
nel T-Ld under different transverse anisotropy is illustrated 
in Fig. 5, and the interatomic interaction is still tuned by 
p-wave Feshbach resonance centered at Bq = 198.8G. 
The vertical dashed lines indicate where the p -wave CIR 
occurs. 


V. CONCLUSIONS 


whose position shows an apparent shift with increasing 
transverse anisotropy. Besides, a two-channel mechanism 
is presented for p -wave confinement-induced resonance 
in a waveguide, which was first proposed for s-wave in¬ 
teraction. We find this effective two-channel picture is 





Figure 5: (Color online) The molecular state of the closed 
channel Hamiltonian Hd for different transverse anisotropy r/ 
. The horizontal dotted lines denote the continuum threshold 
of the open channel, while the vertical red dashed lines in¬ 
dicate the position where a zero-energy scattering resonance 
occurs. 


We theoretically investigate the influence of the trans¬ 
verse anisotropy of the confinement on the two-body 
state with p-wave interaction in a quasi-one-dimensional 
waveguide. The two-body problem in such quasi-one- 
dimensional systems is solved. The interatomic inter¬ 
action is modeled by p -wave pseudopotential. We find 
there are totally three bound states due to the non-zero 
orbital angular momentum of p -wave interaction, while 
there is only one two-body bound state supported by 
s-wave pseudopotential. In addition, the effective one¬ 
dimensional scattering amplitude and scattering length 
are derived. We predict the p-wave confinement-induced 
resonance for any transverse anisotropy of the waveguide, 


still valid for p-wave interaction. All our calculations are 
based on the parameterization of the 40 K experiments, 
and can be confirmed in future experiements. 
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